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The lifetime of ground-state hadrons containing a 6 
quark and lighter quarks is largely determined by the 
charged weak decay of the 6 quark. Interactions involv- 
ing the lighter quarks, referred to as spectator processes, 
alter 6-hadron lifetimes at approximately the 10% level. 
Lifetimes are important to probe our understanding of 
the low-energy strong interaction. While precise predic- 
tions for 6-hadron lifetimes are difficult to calculate, ra- 
tios are predicted with fairly high accuracy by the Heavy 
Quark Expansion (HQE) [T] . This framework of theoret- 
ical calculation is used to predict low energy QCD effects 
in many flavor observables. For example, HQE predicts 
the decay-width of B s mesons to final states common 
to Bg and B®, T s 12 , which enters the decay- width dif- 
ference in the system and several CP violation ef- 
fects. The measurement of lifetime ratios provides a sim- 
ple and accurate way to test the HQE framework as non 
standard model effects are expected to be highly sup- 
pressed in lifetimes. The ratio t(B + ) /t(B°), R + , (charge 
conjugates are implied throughout) is predicted to be 
in the range 1.04-1.08 [IHI]. Predictions for the ratio 
t(A®)/t(B°), i?A, in HQE, which do not presently in- 
corporate next-to-leading order QCD corrections, lie in 
the range 0.88 ± 0.05 [21 HI IS]. The first measurements 
of the A° lifetime have been at the lower end of that 
range. However, recent high precision measurements by 
the CDF experiment [51 [7] based on 1.0 fb _1 of data, 
are significantly higher than previous results. It's there- 
fore useful to keep pursuing lifetime measurements with 
increased precision to settle the issue. In this letter we 
report precise measurements of 6-quark meson lifetimes 
using the channels B + J/tj) K + , B° J/ip K*° , and 
B° -> J/tJj K° , in addition to the lifetime of the Ag 
baryon using the A® — > J/ip A decay channel. Our data 
sample corresponds to an integrated luminosity of 4.3 
fb _1 and consists of pp collisions at a center of mass en- 
ergy Vs = 1.96 TeV collected by the CDF II detector at 
the Fermilab Tevatron. The measurement reported here 
improves the previous CDF measurement [5] of the A° 
lifetime by updating it with significantly more data. In 
all decay modes, the decay position of the b hadron is esti- 
mated using only J/ tp decay products so that differences 
in decay time resolution between channels is reduced and 
certain systematic uncertainties cancel in ratios of life- 
times. 

The components of the CDF II detector relevant to 
this analysis are described briefly here. Charged particles 
are reconstructed using six layers of silicon microstrip 
detectors with radii between 2.4 cm and 23 cm [5] and 
an open-cell drift chamber called the central outer tracker 
(COT) 9J. These are immersed in a 1.4 T solenoidal 
magnetic field and cover the range \rj\ < 1, where 77 is 
the pseudorapidity defined as 77 = — lntan(#/2), and 9 is 
the polar angle [TUj ■ Four layers of planar drift chambers 
(CMU) [TQ detect muons with p T > 1.4 GeV/c within 
\rj\ < 0.6. Additional chambers and scintillators (CMX) 



[12") cover 0.6 < | ^7 1 < 1.0 for muons withpy > 2.0 GeV/c. 

The reconstruction of 6-hadron candidates begins with 
the collection of J/ip — > n + candidates using a dimuon 
trigger. The extremely fast tracker (XFT) [13] uses 
COT hit information to measure the transverse momen- 
tum and azimuthal direction of charged tracks. Events 
with J/-0 — > candidates are recorded for further 

analysis if two or more extrapolated tracks are matched 
to CMU or CMX track segments, opposite-charge and 
opening- angle requirements are met, and the J/ip can- 
didate has mass in the range 2.7 to 4.0 GeV/c 2 . Af- 
ter offline reconstruction, tracks corresponding to two 
triggered muon candidates are constrained to originate 
from a common vertex to make a J/ip — > can- 
didate. To ensure a high-quality vertex for the life- 
time measurement, each muon track is required to have 
at least three hits in the silicon system. The recon- 
structed invariant mass is required to be in the 
range 3.014 < m(fifi) < 3.174 GeV/c 2 . The b hadron 
is assumed to originate from the average beamspot de- 
termined as a function of time using inclusive jet data. 
The primary vertex for a given event is the x — y posi- 
tion of this beamspot at the average z coordinate of the 
muon tracks at their closest approach to the beamlinc. 
The typical beamline size is f=a 30 fim in x — y, and this 
dominates the uncertainty on the decay length. The pro- 
jection of the transverse decay vector onto the 6-hadron 
Pt direction, L xy , and its uncertainty, cr xy , are also ob- 
tained and are used to estimate the proper decay time, 
ct = , and its uncertainty a ct , where M and px are 
the mass and transverse momentum of the b hadron. The 
uncertainties in the primary and J/ip vertices, and the 
transverse momentum are all included in a ct which has a 
typical value around 0.1 ps. Uncertainties in transverse 
momentum have a negligible effect on ct measurement, 
in comparison to the uncertainty on the vertex positions. 

We reconstruct K*° -> K + n~, K° -> tt + it~, and 
A" — > pn candidates from pairs of oppositely-charged 
tracks fit to a common vertex. As and A de- 
cays can occur outside some layers of the silicon sys- 
tem due to their long lifetime, their tracks are not re- 
quired to have silicon hits. The fitted mass is required 
to be in a mass window; for the K*° this window is 
0.84 < m{Kir) < 0.96 GeV/c 2 , (the lower range is se- 
lected in order to avoid reflections from the (j> —> K + K~ , 
where one kaon is misreconstructed as a pion) , for the K® 
it is 0.473 < m(7T7r) < 0.523 GeV/c 2 , and for the A it 
is 1.107 < m(p7r) < 1.125 GeV/c 2 . This corresponds to 
approximately ±3cr, where a is the mass resolution of the 
reconstructed signal. We suppress K® and A cross con- 
tamination by rejecting K® (A ) candidates with proton- 
pion (pion-pion) invariant mass consistent with A (A"°). 
We reconstruct the 6-hadrons by performing a kinematic 
fit of all 6-hadron final state tracks to the appropriate 
topology: two spatially separated vertices in the case of 
A" -> J/iphP and B° -> J/ipK®, one vertex in all other 
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cases. A mass constraint is applied in the J/ip fit, and 
the reconstructed momenta of the K® and A are required 
to point back to the J/ip vertex. We exclude candidates 
with a ct > 100 fim to ensure well measured vertices. Ad- 
ditional selection requirements implying consistency with 
the fit assumptions (common vertex or vertices, mass and 
pointing constraints) are also applied. Further selection 
requirements on the transverse momenta of the 6-hadrons 
and daughter particles, invariant mass of the K° , K*° , 
and A , the vertex probability of the 6-hadrons, and the 
L xy significance of the K® and A were obtained via an 
optimization procedure, which maximizes the quantity 
S/yS + B over all of the selection requirements. The 
number of signal events (S) is estimated from simulation 
and the number of background events (B) from the mass 
sidebands in data. Sidebands are events away from the 
mass peak and form a sample of pure background. 

For B + and B° modes, only candidates with a recon- 
structed B mass between 5.17 and 5.39 GeV/c 2 are used 
for the lifetime measurements. For the A" mode, the 
mass range is set to 5.48 — 5.76 GeV/c 2 . These ranges 
provide a sufficient number of events in the sideband re- 
gions to constrain the background shape while avoiding 
regions where the mass distribution has complex struc- 
ture. The invariant mass distributions for are shown in 
Fig. [I] where the sideband regions are indicated. The 
hadron masses are consistent with world average values. 
We observe the following yields of signal events: 45000 
± 230 (B+), 16860 ± 140 (B° -> J/ip K*° ), 12070 ± 
120 (B° -> J/ipK°), and 1710 ± 50(A°). The life- 
times are extracted using an unbinned maximum likeli- 
hood method. The likelihood function C is multivariate, 
and is based on the probability of observing a candidate i 
with reconstructed mass, m,, decay time, cti, decay time 
uncertainty, of*, and mass uncertainty, cr|™. The PDF is 
factorized in the following form: 

£ = 11^ • P m(miWr) ■ T°(cU\af) ■ S^{af) (1) 

i 

+ (1 - /.) ■ (m,-) • T?(cU\af ) • SU{af )], 

where P m , T ct , and S a ct and are the normalized prob- 
ability density functions (PDF) for observables m i: cti 
and of*, the superscripts s or b refer to the PDF for sig- 
nal or background candidates, respectively, and f s is the 
fraction of signal events. The PDF S^ct is substantially 
different for signal and background events and therefore 
needs to be taken into account as discussed in Ref. [14] . A 
PDF term for cr™ can be ignored as the distribution of a m 
is observed to be similar for both signal and background 
and hence represents a constant in the log-likelihood. 
The signal mass distribution, P^ lll is modeled as 

2 Fisexv ( / _ ( m »- m ') 2 N \ 



where F\ + F2 = 1, mo is the mass of the hadron, and 
Uk are scale factors to account for the misestimation of 
the mass resolutions. We find that two Gaussians are 
sufficient to model the data. The background mass dis- 
tribution, P^, is modeled as a linear function. 

The signal ct distribution, T c s t , is modeled by an expo- 
nential (e _cti / cr /cr) convolved with a detailed detector 
ci-resolution model, TZ. The background ct distribution, 
T c b t , has four components: a S— function convolved with 
TZ to account for backgrounds from prompt J/ip's origi- 
nating from the primary vertex, and one increasing and 
two decreasing exponentials that account for mismea- 
sured decay vertices and background from other heavy- 
flavor decays. These exponential components are con- 
volved with a single Gaussian of width er? multiplied by 
a scale factor. The relative contribution of each back- 
ground component is determined by the data. The pa- 
rameters of the background model are mainly determined 
from the candidates in the mass sidebands. Studies of in- 
clusive 6-hadron decays have shown that after the selec- 
tion requirements the contamination from other b decays 
is very low and, furthermore, that the mass distribution 
of the long lived background components is flat in the 
fitted mass range, and hence the mass sidebands can pro- 
vide a realistic background model for candidates in the 
signal mass range. 

The same resolution model, TZ, is used for signal and 
prompt background events. The detector resolution is 
based upon a Gaussian with width Sj ■ erf*, where Sj 
is the scale factor that accounts for the misestimation 
the of*. Motivated by a study of resolution in an in- 
clusive sample of J/ip events, where prompt J/ip events 
dominate, 1Z is modeled as TZ = Ylj=x hi i^ 7TS j a i t ) ' 
exp (-t 2 /2{ Sj af ) 2 ), where fi + f 2 +fs = 1. The Gaus- 
sians are centered at zero as no evidence of an offset in the 
data samples is observed. Small differences in TZ arise be- 
tween decay channels due to different x 2 distributions for 
the vertex fits of decays with different number of tracks. 
Therefore the parameters fj and Sj are obtained sepa- 
rately for each channel from a fit to data in the mass 
sidebands. This yields an accurate determination of TZ 
since the background events are primarily expected to 
originate from the interaction vertex. 

The functional form of the PDF S a a is determined em- 
pirically using data in the B hadron mass sidebands as in- 
dicated on Fig. [T] The parameters of the function, which 
are different for signal and background are determined 
from the final fit to data. After the resolution model 
parameters are determined from the mass-sideband only 
fit, the likelihood is calculated for each candidate and the 
product is maximized in each of the four channels to ex- 
tract the lifetime, signal yield and the other parameters 
required to describe the mass, background decay time, 
and a ct distributions. Decay time projections of the like- 
lihood function are compared with the data in Fig. [2] 
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FIG. 1: Invariant mass together with mass fit projection for (a) B + -> J/ipK+, (b) B° -> J/ijjK* , and (c) Ag -> J/^A°. 
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FIG. 2: Decay time distributions for (a) B + -> J/ipK + , (b) B° -> J/if)K*, and (c) A° -> J/i/>A° candidates. 



We considered correlated and uncorrelated systematic 
uncertainties. Correlated uncertainties affect all mea- 
sured lifetimes identically, and cancel in ratios. We esti- 
mate uncertainties due to any residual misalignments of 
the silicon detector using Monte Carlo samples generated 
with radial displacements of individual sensors (internal 
alignment) and relative translation and rotation of the 
silicon detector with respect to the COT (global align- 
ment). The XFT triggers on tracks assuming they orig- 
inate from the center of the beam, which may introduce 
a bias for triggering long-lived decays. No indication of 
any bias was found in a study of the XFT response in 
a large sample of simulated events but a small uncer- 
tainty is assigned due to the limited statistical precision 
of the evaluation method. The systematic uncertainty 
that results from ignoring the correlation between re- 
constructed mass and <r ct in the likelihood is found to 
be negligible. The remainder of systematic uncertainties 
arc treated as uncorrelated. They were determined us- 
ing pseudo-experiments in which many statistical trials 
are generated according to alternate PDFs where the al- 
ternate parameters are derived from data. The shift in 
data due to the alternate PDFs were consistent with the 
shift observed with the pseudo-experiments. As the time- 
resolution is determined from the prompt events, and the 
shape of those events is sensitive to the modeling of long- 
lived (positive and negative) background, uncertainties in 
the background modeling can affect the lifetime through 
the resolution function. We account for that uncertainty 
by including an extra long-lived component in the back- 



ground model. This alternate description produces a sub- 
stantial change in the fraction of prompt events (approx- 
imately 7%), and has a small but non-negligible effect 
on the lifetime. A further small uncertainty arising from 
the functional form of 1Z is also assessed and included 
in the total resolution uncertainty. To evaluate uncer- 
tainties in the mass model, alternate parametrizations, 
including a 2nd order polynomial for background, and 
a single Gaussian to describe signal events, were con- 
sidered. Alternatives to the background PDF included 
extra long lived and Gaussian components. We deter- 
mined the uncertainty due to the a ct parametrization by 
using a reasonable alternate model. We also considered 
the effect of ignoring any differences between signal and 
background mass uncertainties by using distributions de- 
termined from data to generate the values of the mass 
uncertainty in the pseudo-experiments. We also deter- 
mined the systematic uncertainty due to the presence of 
the Cabibbo suppressed channel B + — > J/ipir + in the 
charged B decays, and the effect of swapping the kaon 
and pion hypotheses in K*° reconstruction. The possi- 
bility of a systematic biases caused by the a ct and pr 
selection requirements were found to be negligible. The 
systematic uncertainties are summarized in Table |TJ 

We measure t(B+) = [1.639 ± 0.009(stat) ± 
0.009(syst)] ps and t(B°) = [1.507 ± 0.010(stat) ± 
0.008(syst)] ps where the two B° measurements have 
been combined. These results are consistent and of simi- 
lar precision to the leading measurements from Belle [15] 
which are t(B+) = [1.635 ± O.Oll(stat) ± O.Oll(syst)] 
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TABLE I: Summary of systematic uncertainties. 





J/ipK+ (fs) 


J/ip K* u 


(fs) J/ij)K" s (fs) 


J/^A U (fs) 


R+ Ra 


Resolution function 


2.5 


3.5 


3.0 


8.9 


0.0024 0.0061 


Background ct model 


1.0 


2.3 


4.1 


4.6 


0.0017 0.0034 


Mass model 


2.8 


2.8 


2.8 


2.8 


0.0020 0.0017 


Proper decay time uncertainty 


1.7 


1.7 


1.7 


4.3 


0.0010 0.0029 


Mass uncertainty 


3.0 


3.0 


3.0 


3.0 


0.0020 0.0012 


Total uncorrelated 


-1-^9 


zn u. z 


ZL D . o 


-1-117 


n nn/19 n nn7Q 


Alignment 


6.7 


6.7 


6.7 


6.7 




Cabibbo suppressed mode in B + 


0.7 








0.0004 - 


Swapped track assignment in B° 




0.7 








Possible trigger bias 


1.7 


1.7 


1.7 


1.7 




a ct -m correlation 


0.7 


0.7 


0.7 


0.7 




Total 


±8.7 


±9.3 


±9.7 


± 13.7 


0.0043 0.0079 



ps and t(B°) = [1.534 ± 0.008(stat) ± 0.010(syst)] 
ps. The similarities between the decay channels al- 
low for the accurate determination of the ratio i?+ = 
[1.088±0.009 (stat)±0.004 (syst)] which favors a slightly 
higher value than the current average of 1.071±0.009 [2]. 
These results are consistent with the current HQE pre- 
dictions, giving further confidence in this theoretical 
framework, and also provide an accurate test for fu- 
ture lattice QCD calculations. For the A° we measure 
t(A°) = [1.537 ± (stat)0.045± (syst)0.014] ps and R A = 
[1.020 ± 0.030(stat) ± 0.008(syst)]. This measurement is 
the most precise measurement of t(A®) and is consis- 
tent with the previous CDF measurement in this decay 
channel of r(A°) = [1.593i^|(stat)±0.033(syst)] ps [6 
but is more than 2a larger than the world average of 
1. 383^0^48 ps and the previous CDF measurement [7], 
performed on a different decay channel (Afir T ): [1.401± 
0.046(stat)±0.035(syst)] ps. The ratio is also higher than 
the predicted values of 0.88±0.95. In summary, we re- 
port the most precise determination of t(B + )/t(B°). It 
is consistent with other measurements and the predicted 
value which gives confidence in the HQE framework for 
flavor observables. We also report the most precise mea- 
surement of t(A°), which supports a higher value than 
the world average and theory predictions. 
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